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Abstract —This paper investigates the performance of relay 
networks in the presence of hybrid automatic repeat request 
(ARQ) feedback and adaptive power allocation. The thronghput 
and the outage probability of different hybrid ARQ protocols are 
studied for independent and spatially-correlated fading channels. 
The results are obtained for the cases where there is a sum power 
constraint on the source and the relay or when each of the source 
and the relay are power-limited individually. With adaptive power 
allocation, the results demonstrate the efficiency of relay-ARQ 
techniques in different conditions. 

I. Introduction 

Relay-assisted communication is one of the promising tech¬ 
niques that have been proposed for the wireless networks. 
The main idea of a relay network is to improve the data 
transmission efficiency by implementation of intermediate 
relay nodes which support the data transmission from a source 
to a destination. The relay networks have been adopted in the 
3GPP long-term evolution advanced (LTE-A) standardization 
m and are expected to be one of the core technologies for 
the next generation cellular systems. 

From another perspective, hybrid automatic repeat request 
(ARQ) is a well-established approach for wireless networks 
lu-ini. The ARQ systems can be viewed as channels 
with sequential feedback where, utilizing both forward error 
correction and error detection, the system performance is 
improved by retransmitting data that has experienced bad 
channel conditions. Thus, the combination of relay and ARQ 
improves the performance of wireless systems, because the 
ARQ makes it possible to use the relay only when it is needed. 

Due to the fast growth of wireless networks and data- 
intensive applications in smart phones, green communication 
via improving the power efficiency is becoming increasingly 
important for wireless communication. The network data vol¬ 
ume is expected to increase by a factor of 2 every year, 
associated with 16 — 20% increase of energy consumption, 
which contributes about 2% of global CO 2 emissions lfT4l . 
Hence, from an environmental point of view, minimizing the 
power consumption is a very important design consideration, 
and green data transmission schemes must be taken into 
account for the wireless networks CSl-iEol. Moreover, as 
most wireless devices operate with limited battery power, it is 
very important to find ways of maximizing the device lifetime 
by efficiently utilizing the limited power. These are the main 
motivations for this paper, in which we analyze the power- 
limited performance of the relay-ARQ setups. 
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The basic principles of different ARQ protocols are derived 
in Ei-ii. Power allocation in ARQ-based single-user (with¬ 
out relay) networks is addressed by, e.g., ll9l- lfT3l . Also, ETl - 
ll28l study the problem in relay networks. There are a number 
of papers dealing with energy efficiency and power allocation 
in relay-ARQ setups. These works can be divided into two 
categories, as stated in the following. 

In J29l-ll4n. the source and the relay use, e.g., space-time 
codes (STCs) to make a distributed cooperative antenna and 
retransmit the data simultaneously in rounds when the relay is 
active; With an outage probability constraint, im, m (resp. 
El) study the energy efficiency (resp. long-term average 
transmission rate) of STC-based relay-ARQ systems. The 
energy and spectrum efficiency of the basic and hybrid relay- 
ARQ networks are verified in ll^ - ll34l as well. Also, lITSl 
designs a multi-relay-ARQ network using Alamouti codes. 
Assuming the source and the relay to be close, l36l, El 
investigate the throughput of relay networks using different 
ARQ protocols. Optimizing the delay-limited throughput and 
deriving a closed-form expression for the average power 
of the source are addressed by Il38l and |[39l, respectively. 
Considering the incremental redundancy (INR) protocol, m 
studies the performance of the relay-ARQ setups in fast¬ 
fading conditions. Finally, the results of ll40l are extended in 
El, where the system performance is compared with cases 
having only one of the source or the relay active in the 
retransmissions. Implementation of STCs in these works is 
based on the assumption that there is perfect synchronization 
between the source and the relay. 

In Il42l - ll5^ . only one terminal (either the source or the 
relay) is active in the retransmission rounds, as opposed to 
m-M- For instance, m studies the outage-limited energy 
minimization in single-user and relay-ARQ networks. Oppor¬ 
tunistic relaying, rate adaptation and analyzing the energy- 
delay tradeoff curve are considered by El, El and El, 
||46l, respectively, where the direct source-destination link is 
ignored. Also, the throughput, the packet error rate and the 
effective capacity of different ARQ-assisted relay networks are 
studied in Il47l - ll49l . respectively. Power scaling in MIMO and 
cognitive radio relay-ARQ networks is addressed in ll50l and 
El, respectively. Finally, ll5^ studies a relay-ARQ network 
using superposition coding. References EH-Ell, El, El, 
El-lESl are based on the assumption that there is a fixed 
transmission power for the source and the relay. Meanwhile 
El, El, El, El optimize the power allocation between 
the source and the relay under a sum power constraint, while 
they use the same powers in all retransmissions. Also, ll42l 
investigates the power allocation between the retransmissions 
for basic ARQ schemes and ll^ studies the average power 
of the source with repetition time diversity (RTD) ARQ and 
a fixed power for the relay. 
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In theoretical investigations, the communication links be¬ 
tween the source, the relay and the destination are normally 
assumed to be independent 1^ - 15^ . This is an appropriate 
model for many practical scenarios and makes it 

possible to analyze the system performance analytically. How¬ 
ever, the independent fading channel is not always a realistic 
model. For instance, the relay is normally located close to the 
destination in moving-relay systems m, 1551. As a result, 
there might be considerable correlation between the source- 
relay and the source-destination fading coefficients. Also, e.g., 
l52l demonstrates the cases where the source is connected to 
the destination through a relay which is close to the source. In 
this case, the source-destination and the relay-destination links 
may be spatially-correlated. For these reasons, it is interesting 
to extend the independent fading model to the case where there 
is spatial correlation between the channels. 

In this paper, we study the throughput and the outage 
probability of the relay-ARQ networks in cases where there 
is either a long-run sum power constraint on the source and 
the relay or when each of the source and the relay are power- 
limited individually. Adaptive power allocation between the 
retransmissions is used to improve the system performance. 
We derive closed-form expressions for the average power, the 
throughput and the outage probability of different relay-ARQ 
protocols in the cases with independent or spatially-correlated 
fading channels. Moreover, we investigate the effect of fading 
temporal variations on the data transmission efficiency of the 
relay-ARQ systems. 

As opposed to M-M, we study the scenario where only 
one of the source or the relay is active in each ARQ-based 
retransmission round. Also, the problem setup of the paper is 
different from the ones in Ell-Ell because 1) we consider 
adaptive power allocation between retransmissions of hybrid 
ARQ protocols, 2) the results are obtained with different 
sum and individual power constraints on the source and the 
relay and 3) we investigate the system performance in both 
independent and spatially-correlated fading conditions, with 
noisy/noise-free feedback signals. Finally, our discussions on 
the users’ message decoding probabilities (Theorems 1-3) have 
not been presented before. 

The results show that there is a structural procedure to 
study different performance metrics of relay-ARQ networks 
experiencing different fading models. Optimal power alloca¬ 
tion is shown to be very useful in terms of outage probability, 
throughput and coverage region of the relay-ARQ network, 
when there is a sum power constraint on the source and 
the relay. With individual power constraints on the source 
and the relay, however, optimal power allocation increases 
the throughput (resp. reduces the outage probability) only 
at low (resp. high) signal-to-noise ratios (SNRs). Compared 
to the fixed-length coding scheme, the throughput of the 
relay-ARQ network increases when variable-length coding is 
utilized. With the practical range of spatial correlations, the 
performance of the relay-ARQ network is not sensitive to the 
spatial correlation. However, the data transmission efficiency 
of the network is reduced at highly-correlated conditions. 


II. System model 

We consider a relay-assisted communication setup con¬ 
sisting of a source, a relay and a destination. The channel 
coefficients in the source-relay, the source-destination and 
the relay-destination links are denoted by h®”, and 
respectively. Also, we define p®” = and 

which are referred to as the channel gains in 
the following. A maximum number of M ARQ-based retrans¬ 
mission rounds is considered, i.e., the data is (re)transmitted 
a maximum of M -f 1 times. Moreover, we define a packet 
as the transmission of a codeword along with all its possible 
retransmission rounds. In each packet, Q information nats are 
sent to the destination and the length of the subcodeword used 
in the m-th round of the ARQ is denoted by Im- Thus, the 
equivalent data rate, i.e., the code rate of the ARQ, at the end 
of the m-th round is given by i?(m) = ; . 

We study the system performance for two different block¬ 
fading conditions: 

• Quasi-static. In this model, the channel coefficients are 
assumed to remain fixed within a packet period, and then 
change to other values based on their probability density 
functions (pdf). 

• Fast-fading. Here, the channel coefficients are supposed 
to change in each retransmission round. 

The quasi-static model, studied in Subsections IV.A-B, 
represents the scenarios with slow-moving or stationary users, 
e.g., ED, im, Ea. On the other hand, the fast-fading, studied 
in Subsection IV.C, is an appropriate model for the high speed 
users and frequency-hopping setups where the channel quality 
changes in the retransmissions independently, e.g., ED, ED, 

Ei- 

In each link, the channel coefficient is assumed to be known 
by the receiver, which is an acceptable assumption in block¬ 
fading channels ||9ll- E3l . However, there is no instantaneous 
channel state information available at the transmitters except 
the ARQ feedback bits. The ARQ feedback signals are initially 
assumed to be received error-free, but we later investigate the 
effect of erroneous feedback bits as well (Section IV). 

Relay-ARQ model: The considered relay-ARQ protocol 
works as follows. In each packet period, the data transmission 
starts from the source. If the data is decoded by the destination, 
an acknowledgement (ACK) is fed back by the destination to 
the source and the relay, and the retransmissions stop. Oth¬ 
erwise, the destination transmits a negative-acknowledgment 
(NACK). Only one terminal (either the source or the relay) is 
active in each retransmission round; the relay becomes active 
and the source turns off, as soon as the data is decoded by the 
relay. That is, if the relay successfully decodes the message, it 
sends an ACK to the source and starts retransmission until 
the destination decodes the data or the maximum number 
of retransmissions is reached. In other words, with error-free 
feedback bits, the following cases may occur during a packet 
transmission period: 1) receiving an ACK from the relay and 
a NACK from the destination, the source turns off and the 
relay starts retransmission. 2) With NACKs from the relay 
and the destination, the data is retransmitted by the source. 
3) Receiving an ACK from the destination, the source ignores 
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the ARQ feedback of the relay and the retransmissions stop 
(Performance analysis in the cases with noisy feedback bits is 
studied in Section IV.D.)- 

The motivations for considering the proposed data transmis¬ 
sion model are as follows. Letting the relay retransmit instead 
of the source when the source-destination link experiences bad 
condition makes it possible to exploit the potential diversity 
gain through the relay channel. Also, in practice, the relay 
is located such that the relay-destination link experiences 
better average characteristics than the source-destination link. 
Therefore, it is more beneficial to use the power resources for 
the relay, instead of dividing the power between the source and 
the relay, if the relay decodes the message. Finally, as seen 
in the following, for Rayleigh-fading conditions the proposed 
scheme outperforms the state-of-the-art approaches, in terms 
of outage probability/throughput. 

III. Problem formulation 


Definitions 

r-T 


Deriving the_^expected values 


Calculating the probabilities 
Pr(AJ, Pr(B„.,), Pr(S J 


Closed-form expression for 
the functions of interest 


Optimization algorithms 



Simulation results 


Parameter optimization and 
simulations 


Figure 1. An overview of the paper. First, we obtain closed-fonn expressions 
for the functions involved in (T). Then, iterative optimization algorithms are 
applied to optimize the parameters based on the closed-form expressions. 


In this paper, we study the problem of 
max O, O = Ip, — Pr(Outage)| 

subjectto A, A = ($’* < 0®)&((I>‘' < cj/)}. 

( 1 ) 

In words, we investigate the long-term throughput p and the 
outage probability Pr(Outage) as the evaluation yardsticks. 
The optimization parameters are the equivalent data transmis¬ 
sion rates i?(m) as well as and m = 1,M + 1, 
which denote the source and the relay power used in the 
m-th retransmission round, respectively (because the noise 
variances are set to 1, P^ and P^, in dB 101og2^Q(P^) and 
10log]^Q(P,!i^), represent the transmission SNR as well). Fi¬ 
nally, the throughput and the outage probability are optimized 
under two different power-limited scenarios: 

• Scenario 1. The total power for data transmission in 

the relay-ARQ setup is limited, which is represented by 
(jjtotai £ (^totai (j)totai jj^g power in the source 

and the relay, averaged over many packet transmissions, 
and denotes the total power constraint. This scenario 
is of interest in the green communication concept, where 
the goal is to minimize the total average power required 
for data transmission QSl-llol, and also for electricity- 
bill minimization. 

• Scenario 2. There are individual power constraints on 
the source and the relay, which is represented by (<I>® < 

< (j/) in ([T]i. Here, and are the average 
power in the source and the relay, respectively, and </>* and 
(jf denote their corresponding thresholds. This scenario 
models the case where the source and the relay are 
battery-limited 13-113. 

To study ([T]i, the following procedure is considered (please see 
Fig. 1 as well). First, we derive closed-form expressions for the 
functions rj, Pr(Outage), and which are involved 

in O- Then, since O is a nonconvex problem, iterative 
optimization algorithms are used to optimize the parameters 
based on the closed-form expressions. 


In three steps we obtain the closed-form expressions of 
rj, Pr(Outage), and The first step is to de¬ 

fine the metrics and the constraints as functions of a few 
expected values. Then, in the second step, we derive the 
expectations as functions of predefined probability terms. 
The last step is to represent the considered probabilities 
in terms of P^, P^, m = -f 1, i.e., the 

optimization parameters of O- Interestingly, the two first 
steps are independent of the considered ARQ protocol and 
the fading channel model. Thus, they are explained in the two 
following subsections. The third step, however, depends on the 
characteristics of the ARQ schemes and the fading channel 
model. For this reason, we specify the results for different 
ARQ protocols and fading channel models in Sections IV and 
V. 


A. Step 1: Definitions 

The outage probability is defined as the probability of the 
event that the data can not be decoded by the destination when 
the data (re)transmission is stopped. Also, the throughput (in 
nats per channel use (npcu)) is given by a, El, ini, ina 

Y.k=iQk 

n = lim —- 

K-^oo 


= lim 
K—^oo 


K^k=l Qk ig) E{Q} 


( 2 ) 


Here, Qk is the number of information nats successfully 
decoded by the destination in the fc-th packet transmission. 
Also, denotes the total number of channel uses in the 
fc-th packet transmission, i.e., if the message 

retransmission of the fc-th packet continues for m rounds 
(see ©). Note that in each packet part of the data may be 
(re)transmitted by the source or the relay and 
where t\. and are the source and the relay activation periods 
in the fc-th packet transmission, respectively. In general, 
and are random values which follow the random variables 
Q and respectively, as functions of the channel realiza¬ 

tions. Also, (a) in (|2|l is based on the law of large numbers. 
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with E{.} representing the expectation operator, and the fact 
that the limits, e.g., lim ^ YlE Qk, j™ Ef=i 
exist Q, HU, m. 

With the same arguments, the average power terms 
<!)'* and are obtained by 


= lim 
K—¥00 


= lim 
K—^oo 


E 

E 

E 

E 


K 

k^l 

K 

k^l 


K 

k^l 

K 

/c=l 


a 

a 


Em 

E{r^}' 


Em 

E{V}^ 


(3) 

(4) 


Ctotal I cr 

_ _ y Z^fc=l >fc Z^k—1 Sfc 

^ K^oo 

l^km ^k l^km ^k 

Em+Em ... 

E{rm ■ ' 

Here, Ck are the source, the relay and the 

total transmission energy in the fc-th packet transmission, 
respectively, with = Ck + Ck- Also, S'", T'* and T‘ 
denote the random variables corresponding to and 

respectively. Note that the metrics and constraints are 
functions of a few expected values. 


If the data is decoded at the end of the m-th round, the total 
number of channel uses is l(^m) = EEi Also, the total 
number of channel uses is Z(m+i) = E^T^ if ™ outage 
occurs, where all possible retransmission rounds are used. 
Thus, the expected number of total channel uses, i.e., 
in (| 2 ]i and (|5]l, is obtained by 


M+1 / m \ 
m—1 \i—l / 


M+1 

Pr(Am) + ( ^ k) Pr(Outage). 
mi 


(8) 


From (|7]i-(|8]l and i?(m) = ^+ 7 -, the throughput (|2]) is found 
as 

EEtPr(A^) . 9 . 

' y^M+l Pr(A,n) l-E+IPrC^m) ^ ' 

2—im=l R^rn) fl(M + l) 

If the source stops data (re)transmission at the end of 
the m-th round, the total energy consumed by the source is 
^\m) ~ EEi Elh- Therefore, the source consumed energy is 
a random variable given by 

m 

= if m = l,...,M + l, ( 10 ) 

i=l 


B. Step 2: Deriving the Expected Values 

Let us define the following events: 

• Am is the event that the data is successfully decoded by 
the destination in the m-th (re)transmission round while 
it was not decodable before. In this case, the codeword 
may have been sent by the source or relay. 

• Pn,m represents the event that the relay is active in rounds 
n + 1,..., m with n = 1,..., M, n < m. In this case, 
the source message has been decoded by the relay in the 
n-th round and, consequently, the source turns off in the 
successive retransmissions. The relay data retransmission 
is stopped in the m-th round if the destination can decode 
the data or the maximum number of retransmissions is 
reached. 

• Sm is the event that the source stops data retransmission 
in the m-th round. In this case, either the maximum 
number of retransmissions is reached or the data has been 
decoded by the relay or the destination. 

The defined events are used to express ©-(IS]) as functions 
of Pm Pirn The details are 

explained as follows. 

According to the definitions, the outage probability is found 
as 

M+1 

Pr(Outage) = 1 — Pr(Am). (6) 

m—1 

If the data is decoded by the destination at any (re)transmission 
round, all Q information nats of the packet are received. 
Hence, the expected number of received information nats in 
each packet is 

M+1 

E{Q} = Q{1- Pr(Outage)) = Q Y. (7) 

m—1 


and, using li = ~ 


M+1 / / m 


Em = Y 


m—1 \ \ 2=1 
M+1 / / m 




m—1 \ \ 2=1 


^{i) -^(z — 1 ) 


Pr(5™) j . 
( 11 ) 


With the same arguments, the expected activation period of 
the source, i.e., E{T^} in Q, is found as 


M+1 


M+1 


E{r} = E E^O = Q E 


Pr(5™) 


m—1 \ \i—l 


m—1 


E(^rn) 


which, along with (fTTIi . leads to 

„ _ E";; ((E”. p/ (nb - %b7)) 


= 


E 


M +1 Pr(5m) 
R(m) 


( 12 ) 


(13) 


Given that the data is retransmitted by the relay in the 
n + l,...,m rounds, its consumed energy is = 

EEn+i Elk which is consumed during E™ n+i k channel 
uses. Thus, we can use the definition of Bn^m, the same 
procedure as in ([T0ll-([T3ll and J2jLn+i k = - -r^) 

to write 


E{^ } — C(n+l,m) Pl'( + ,m) 

Vn<m,m< M+1 


= <3 E 

Vn<m,m< M+1 



1 



( 14 ) 
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m 

E{r}= Y. (E 

Vn<m,m<M+l 2 =n+l 

= Q 'Y/ 

Vn<m,m<M+l 



= 


m'} 

E{V} 


^Vn<m,m<M+l n +1 ^■;( 7 fe- 7 iT 7 ))P'(«-)) 


E 


Vn<m,m<M+l 




(16) 


Note that the summations in dH-lUlll are on all possible 
activation conditions of the relay. Finally, from (|8]l, (fTTl i. 
(fT4li . the total transmission power ^ is 

obtained by 


^total _ _ _ _ 

“ y^M+l Pv{Am) , 1-E^=1 Pl-(^m) ’ 

2^m=l Ji(M+l) 



+ E 

Vn<m,m<M+l 



PriSm) 




Pr(S„,m) 

(17) 


From it follows that the only difference be¬ 

tween different ARQ protocols is in the probability terms 
Pr(Am), Pr(S'm) and Pr(i7„_m)- Also, to derive closed-form 
expressions for the outage probability, the throughput and 
the average power terms, the final step is to represent the 
probabilities Pr(Am), Pr(5'm) and Pr(i?„„j) as functions 
of -Pm I Pm: m = i.e., the optimization 

parameters of O- Sections IV and V are devoted to obtain 
the probability terms for different ARQ protocols and fading 
channel models. 


IV. Performance analysis in 

SPATIALLY-INDEPENDENT RAYLEIGH-FADING CONDITIONS 

For the spatially-independent Rayleigh-fading channels the 
fading coefficients follow ^ CAf(0, ^), ^ CM(0, 

and ~ CN{Q, Thus, the pdf of the channel gains 
are given by /g:r(a;) = fg«i{x) = and 

fgr<i{x) = Here, A*"^, A“^ and A***^ are the fad¬ 

ing parameters determined by the path loss and shadowing 
between the corresponding terminals. Performance analysis 
of the relay-ARQ setup in the presence of RTD and INR 
protocols is as follows. 


the initial data rate and the length of the codeword, 
respectively. Moreover, with MRC, the received SNR of, 
e.g., the relay at the end of the m-th retransmission round 
increases to EEi Thus, the data is decoded by 
the relay at the end of the m-th round (and not before) if 
log(l + 5 - j:TY P!)<R< log(i + YZi Pl)^ which is 
based on the fact that with SNR x the maximum decodable 
rate is ^ log(l + x), if a codeword is repeated m times ||7l . 
ifT^ . In this way, Pr(5m), i-C-, the probability that the source 
stops retransmission at round m, is found as 


Pr(5m) = 


dm “t“ /^m if ur — 1, . . . , ilT, 
7m if m = M -I- 1, 

m —1 

= Pr(log(l + p-E^')<^ 


< log(l + Y ^ E 

m—1 

Pm = Pr(log(l + Yp')<P 

m m— 1 

< log(l + 5^'* E ^ E 


2=1 

M 


M 


7m 


= Pr(log(l + E ^ ^ E < py 


(18) 


Here, ctm is the probability that the relay decodes the data at 
round m, before the destination. Moreover, /3m denotes the 
probability that the destination decodes the data at round m 
while the relay had not decoded the message up to the end 
of the (m — l)-th round (the message may be decodable by 
the relay at the m-th round). Also, the source retransmits the 
codeword M+1 times, if none of the relay and the destination 
have decoded the data until the M -th round, leading to jm in 
(fTSl) . Note that J2m=i Pr(S'm) = 1, as a maximum of M -f 
1 (re)transmissions is considered. For independent Rayleigh- 
fading channels, (OS is rephrased as 




FA 




9' '•■^m-l ps 
Z^i=l 


- 1 f 

:)-Fg4 


- 1 


E m p; 
2=1 


7 ) FA 


- 1 


E ’TI ps 

2=1 2 




), 


Pm 



)-Fg4 


- 1 


E m ps 
2=1 


) Fg.{ 




- 1 


Y^m-l ps 

Z^i=l 


:) 



A. RTD Protocol in Quasi-Static Conditions 

Using the RTD protocol, the same codeword is 
(re)transmitted in each round and the receiver performs 
maximum ratio combining (MRC) of the received signals. 
Thus, the equivalent data rate at the end of the m-th 
round is i?(m) — '§d — m with R and I representing 



With the same procedure, the probability that the destination 
decodes the codeword at the m-th round (and not before), i.e., 
Pr(Am), is obtained by 
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J = 1 


J -1 


= Pr ( log(l + 5 ^^ ^ < i? < log(l + g- ^ PI) n 

i=l i=l 

0 m-l 

\og{l+g^^Y.P^i+9^" Y. PI) 


i=j+l 


< log( 1 + 5^“* Yp^ + 9'“ E ^>^)) ■ (20) 

i—1 i=j-\-l ^ 


Here, „j is the probability that the relay decodes the 
codeword at the j-th round and helps the destination until it 
decodes the message at round m, j < m. Thus, ( l20b gives the 
message decoding probability of destination for all possible 
activation conditions of the relay. For independent Rayleigh- 
fading condition, Ej^m is found as 


active in rounds n + 1 ,..., m, is determined as 


R i if TO ^ M + 1, 

- I if m = M + 1, 

y n— 1 n 

= Pr ( log(l + Yp^)<P^ + 9" E P^) 


2=1 


2=1 


M 


nlog(l + /‘'^P? + 5'‘' ^ pI)<r) = 


Pn — 


e-R-l 
■“eJLi Pi 


M 


2=1 k—n-\-l 

' .R 


^nPrii 


_/^sd (^x')Fgrd 


pr 

A^k—n+1 ^ k / 


dx 


= Qn{ E 


( 22 ) 


2 =n+l 


where £n,m and Qn{x) are defined in (l20l i and (EB, respec¬ 
tively, and Pn is obtained with the same procedure as in (|2TI) . 

Using (fT8])-(l22l). we can express the outage probability, the 
throughput and the average power functions of the relay-RTD 
scheme in terms of R(^rn), Pm ^ Prm to = and 

investigate the system performance, as stated in the following. 


e«-l 


= R ( e“-l t _ R e -1 

\ ^j-1 ps / ^E'^ 

‘ _Y- 


B. INR Protocol in Quasi-Static Conditions 


-A'" 


= e 


— e 




= /o 


/ sd(a;)x 


Considering a maximum of M -b 1 INR-based retransmis¬ 
sion rounds, Q information nats is encoded into a parent 
codeword of length /(m-i-i) = 'YYi Then, the code¬ 
word is punctured into M -b 1 subcodewords of lengths 
Pj. (gti pr < g-R _ X — P^x < g^'^ TO = 1 ,..., M -b 1, which are sent by the source/relay 

^ * » 1 * * i-l-i successive retransmission rounds. In each round, all 

received subcodewords are combined by the receivers (relay 
and destination), to decode the message. In this case, the 
results of ll57ll . Il58l chapter 15], ||59] chapter 7] can be used 
to show that the maximum data rates which are decodable by 
the relay and the destination at the ru-th round are obtained 
by 
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( 21 ) 


where ujj is the decoding probability of the relay at round j 
(and not before). Also, 9j^m represents the decoding probabil¬ 
ity of the destination at the ru-th round, given that the relay 
is active in rounds j -b 1 ,..., to. 

Finally, Pr(i3„ m), i-e., the probability that the relay is 


uUig^^gZ 

ELi h iog(i + A^P!) + EZj+i h iog(i + A'^Pl) 
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respectively, where (l24l i is based on the assumption that the 
relay is active in rounds j + 1,... ,m. Also, 


YY 


i\9 ) 7 


= R( 




R{i) R{i-i) 


)log(l + /V?) 


(25) 


denotes the maximum decodable rate of the destination at the 
m-th round, given that the relay is inactive. 

Although having high throughput and low outage proba¬ 
bility, variable-length coding INR results in high pocketing 
complexity ina, 113. In order to reduce the complexity, 
fixed-length coding INR scheme is normally considered where 
setting Ijn = in (I25ll-(l25]l leads to R^m) = 

^ m 

f;^fixed-le„gth(^sr) ^ ^ ^ (26) 

m 


^d, fixed-length 




( sd id\ 
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2=1 
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=i-ti 


= — > iog(i + g" P!)+ >^ iog(i + / PD ,j <m. 


(27) 

From (|23i-(|27li, we can find the probabilities 
Pr(Am), Pr{Bn,m), Pr(<S'm) for the INR protocol; Replacing 
the terms, e.g., log(l -I- of ’^he RTD by 

R{m) E™ 1 (- R^) log(l + g'^'Pi) for the INR, one 
can use the same procedure as in (fT8l l- (l22l i to recalculate the 
parameters Pm, 7 m , Wj, f^j,m, Pn and, consequently, 
Pr(Am), Pr(i?„ m), Pr(S'm) for the INR. The details 
are presented in the Appendix where the probabilities 
are determined by obtaining Pr(t/^(p'*'') < R(m)) and 

As explained in the Appendix, two-dimensional numerical 
integrations should be used to determine the probability terms 
Pn and 9j^m for the INR, which are difficult to find. This 
is particularly because the boundaries of the two-dimensional 
integrals can not be expressed in closed-form. For this reason, 
upper bounds of and Oj m are presented in the Appendix 
which are tight at low SNRs. Moreover, Theorems 1-2 provide 
other approximation methods which simplify the performance 
analysis in the presence of the INR protocol. The good 
point in Theorems 1-2 is that the two-dimensional integrals 
with unknown integration boundaries are replaced by either 
closed-form expressions or one-dimensional integrals having 
known boundaries. As a result, the terms dj^m and can be 
determined easily. 

Theorem 1: For the fixed-length INR protocol, the perfor¬ 
mance of the relay-ARQ setup is underestimated, i.e., the 
throughput is lower bounded and the outage probability is 
upper bounded, via the following inequalities 


Pr u: 


TV, fixed-length 


VTTE^ 




= 1 — e 




(28) 




<1-14 


-S. 1 

gm — 1 




(29) 


where Vj^rn{v) = /p°° ^"da;. 

Proof: Please see the Appendix. ■ 

Due to properties of the Minkowski’s inequality, the bounds 
are tight at high SNRs. Finally, we close the discussions with 
the following theorem which provides an upper-estimate of 
the system performance. 

Theorem 2: For sufficiently low SNRs, the performance of 
the relay-INR protocol is upper-estimated via the following 
inequalities 
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(30) 


P's/ / 


uua+ipmnZj+A^ + iPD^) 


(31) 


Proof: Please see the Appendix. ■ 

The tightness of the bounds is verified in Subsection IV.D. 


C. Performance Analysis in Fast-Fading Conditions 

Considering the fast-fading modelfl e.g., lfT3 . BTI . Il5^ . 
the outage probability, the throughput and the average power 
functions, i.e., (Ell-®, are obtained with the same procedure 
as before while the achievable rate terms, e.g., (|27] | for the 
INR, are replaced by 

TT-d, fast-fading 

^ 3,m 

1 / 7 m \ 

= “ E + sfPi) + E + sfPD ,j <m. 

^ \i=l i=i-l-l / 

(32) 

Here, gf^ and gl’^ are the source-destination and the relay- 
destination channel realizations at the i-th round. Changing 
the achievable rate terms, and recalculating the probabilities, 
is the only modification required for the fast-fading condition 
and the rest of the procedure does not need to be changed. 
Specifically, Theorem 3 provides a method for calculating the 
probabilities in the fast-fading conditions. 

'Under fast-fading channel conditions, the INR protocol is studied with 
fixed-length coding because the length of the codewords is the same as the 
fading block length. 





















Theorem 3: The power-limited throughput/outage probabil¬ 
ity of a relay-ARQ setup utilizing RTD and INR protocols can 
be analyzed via the equalities 
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(34) 


i=i-|-l 


respectively, if the channel is fast-fading. Here, TL^^n [•] 
denotes the generalized upper incomplete Fox’H function 1^ . 
Proof: Please see the Appendix. ■ 

Finally, to enjoy the practical benefits of the relay-ARQ the 
channel code should satisfy the following requirements; 

1) A parent code that can be punctured into rate-optimized 
subcodewords and 

2) decoders at the relay/destination with performance close 
to (fT^ - (l25] l. (l32] | for all retransmissions. 

There exist several practical code designs, e.g., ||2l, S, that 
satisfy these requirements. Moreover, to implement adaptive 
power allocation, the source and the relay should be equipped 
with adaptive power amplifiers. However, as the power adap¬ 
tation is based on the long-term channel statistics with finite 
levels of transmission powers, the power amplifiers can be 
efficiently designed. 


D. Simulations and Discussions 

Using the same arguments as in ifTTIl . ifT^ . ibH . it can be 
showed that both the power-limited throughput maximization 
and the outage probability minimization of ARQ protocols 
are nonconvex optimization problems, even for the single- 
user (without relay) setup. Therefore, the problems should be 
solved via iterative optimization algorithms. In our simula¬ 
tions, the number of optimization parameters is low enough 
to use exhaustive search, which is what we have used for our 



Source input SNR 101ogj|j((|)b (dB) 


Figure 2. Throughput vs the average transmission SNR of the source 
101ogjQ(())®). The results are obtained with individual power constraints on 
the source and the relay (scenario 2) with d?' = lOdB. RTD ARQ protocol, 
A“ = = 0.5, = 1. 


simulations. Also, for faster convergence, we have repeated 
the simulations by using the iterative algorithm of lfT2ll . and 
by using “fminsearch” and “fmincon” functions of MATLAB. 
Using the closed-form expressions, the results have been 
obtained for different initial settings and we have tested the 
fmincon function for “interior-pointf “active-set” and “trust- 
region-reflective” options of the optimization algorithm. In all 
cases, the results are the same with high accuracy, which is an 
indication of a reliable result. Finally, note that, except for Fig. 
6 where we verify the tightness of the bounds in Theorems 1- 
2, the results of the figures are obtained both analytically and 
via Monte Carlo simulations which lead to the same results. 
Thus, Figs. 2-5, 7-15 represent both the analytical and the 
Monte Carlo-based simulation results. 

In the simulations, fixed-length coding is considered for the 
INR protocol, unless otherwise stated. Also, in all figures, 
except Figs. 12-13, the results are obtained for the quasi-static 
conditions. The fast-fading case is considered in Figs. 12-13. 

Performance analysis with individual power constraints 
on the source and the relay: In Figs. 2-3, we study the 
system performance in scenario 2, where there are individual 
power constraints on the source and the relay. The results 
show that 1) with individual power constraints, optimal power 
allocation increases the throughput at low SNRs. However, 
the gain of optimal power allocation is negligible, when the 
goal is to maximize the throughput at high SNRs (Fig. 2). 
2) Considerable outage probability reduction is achieved by 
optimal power allocation, particularly at high SNRs (Fig. 3). 
Finally, 3) increasing the maximum number of retransmissions 
leads to marginal throughput increment, especially at high 
SNRs, when the source and the relay are individually power- 
limited (Fig. 2). 

Figures 4-15 present the simulation results for the case 
with a sum power constraint on the source and the relay, i.e., 
<j)totai ^ ^totai jjj where 10 log]^o('/’'°'^*) denotes the input 
SNR. 

On the coverage region of the relay-ARQ network: Let us 
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Source input SNR 101og^^((t)^) (dB) 

Figure 3. Outage probability vs the average transmission SNR of the source 
10 logio The results are obtained with individual power constraints on 
the source and the relay (scenaiio 2) with = lOdB. Initial transmission 
rate R = 0.5, M = 1, = 0.5, = 1. 


Optimal power allocation """ INR, relay 



Figure 4. Relative coverage region, compared to the single-user (without 
relay) setup without ARQ. Scenario 1 (sum power constraint in 

o; with 10 being the input SNR). M = 1, initial transmission 

rate R= 1. For relay-assisted channel, A“ = = 0.5. 


define the coverage region as 

7^(e|A'^^ A“', M) = {A^'^IA'^ A“, M, Pr(Outage) < e}. (35) 

That is, (iTSl l dehnes the set of fading coefficients A®‘* (for a 
given SNR, M, A“^, and A“) which leads to Pr(Outage) < e. 
Fig. 4 demonstrates the relative coverage region of different 
communication setups, compared to the single-user setup with¬ 
out ARQ, i.e., LD = ■ In other words, each curve 

in Fig. 4 shows the gain of the considered scheme compared 
to the single-user system without ARQ. The higher the curve 
is, the wider the coverage region il|. The results show that, 
compared to the single-user setup, the implementation of the 
relays leads to considerable coverage region increment. Also, 
compared to uniform power allocation (which corresponds 
to, e.g., 1171-11521), our proposed power-optimized scheme 
increases the coverage region of the relay-assisted networks 
substantially. For instance, consider the INR ARQ, M = 1 
and the coverage threshold e = 10“^. In this case, the 
implementation of power-optimized ARQ in single-user and 
relay-assisted setups increases the coverage region by 17 and 
27 times, respectively. With uniform power allocation, there is 
a (almost) hxed gap between the performance of the INR and 
RTD protocols. With optimal power allocation, the difference 
between the coverage regions of the RTD- and INR-based 
schemes decreases. 

On the effect of variable-length coding: The effect of 
variable-length coding on the performance of the relay-INR 
approach is investigated in Fig. 5. Compared to hxed-length 
coding, considerable (resp. marginal) throughput increment is 
achieved by variable-length coding at moderate (resp. low) 
SNRs. Also, the effect of variable-length coding increases with 
the maximum number of retransmissions, i.e., M. 

On the tightness of the proposed bounds: In Fig. 6, we 
compare the throughput bounds achieved via Theorems 1-2 

^Note that the variances of the fading coefficients, that can quantify the 
distances between the transmission endpoints, are given by, e.g., 



Input SNR 101ogj^((|)”“') (dB) 

Figure 5. Effect of variable-length coding on the throughput of the relay-INR 
protocol, uniform power allocation, A“ = A'^‘* = 0.5, A*** = 1. 


with the exact throughput obtained through (|26li-(|27li. The 
results are obtained with optimal power allocation and for a 
hxed initial transmission rate R = 1. As shown, the bound 
presented in Theorem 1 (resp. Theorem 2) is tight at high 
(resp. low) SNR and their tightness decreases as the input 
SNR decreases (resp. increases). 

Throughput and outage probability in power-optimized 
relay-ARQ scheme: Considering the RTD, Fig. 7 demonstrates 
the throughput versus the outage probability. Moreover, Fig. 
8 studies the effect of optimal power allocation on the outage 
probability. Considering the hgures, it is deduced that 1) with 
a sum power constraint on the source and the relay, optimal 
power allocation improves both the throughput and the outage 
probability considerably. 2) For a given outage probability, 
increasing the number of retransmissions leads to substantial 
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—Exact value via simulations 
— Lower bound obtained via Theorem 1 



INR ARQ protocol, optimal power 
allocation, M=l, R=1,X^^=I, X^^=X^'^=0.5 
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— Upper bound obtained via Theorem 2 
— Exact value via simulations 



INR ARQ protocol, optimal power 


-3 


Input SNR lOlogj^Cct)^^'^^) (dB) 



Throughput (npcu) 


Figure 6. On the effect of the bounds introduced in Theorems 1-2, optimal 
power allocation, scenario 1, M = 1, initial transmission rate i? = 1, = 

1, A**' = = 0.5. 


Figure 7. Throughput vs the outage probability, RTD ARQ protocol, input 
SNR lOlogio(0‘°‘“') = 5dB, scenario 1, A'*'^ = 1, A"' = A""* = 0.5. For a 
given outage probability, the initial transmission rate (and the retransmission 
powers, if required) are optimized to maximize the throughput. 


throughput increment, if the (re)transmission powers are allo¬ 
cated optimally (Fig. 7). Intuitively, this is because with more 
number of retransmissions the relay gets more involved and 
the good characteristics of the relay-destination link are prop¬ 
erly exploited. 3) The difference between the outage-limited 
throughput of optimal and uniform power allocation schemes 
increases with the number of retransmissions and decreases 
with the outage probability (Fig. 7). 4) The difference between 
the outage probability of the single-user and relay networks, 
i.e., the gain of implementing the relay node, increases with 
the input SNR (Fig. 8). 

Comparison with the state-of-the-art schemes: In Fig. 9, 
we compare the data transmission efficiency of the proposed 
relay-ARQ approach with ones in different data transmission 
approaches. As demonstrated in the figure, the proposed relay- 
ARQ approach outperforms the other schemes including 1) 
the single-user channel without ARQ lfT2l . 2) single-user setup 
using ARQ ifTTI . 3) the relay channel without ARQ, e.g., 121], 
El with one relay, 4) the relay-ARQ in the cases where 
the source-destination link is ignored, e.g., HD, ia with 
one relay, or 5) the relay-ARQ when, using STC, the data is 
simultaneously retransmitted by both the source and the relay 
in the retransmissions Ea, iOl, ED. Moreover, the gain of 
the proposed scheme, compared to the considered state-of-the- 
art protocols, increases with the SNR. 

Here, it is interesting to note that, compared to the single- 
user setups, the implementation complexity increases in relay- 
assisted systems; this is because the data transmission is based 
on more handshakings between the terminals when the relay 
is utilized. Also, compared to the single-user systems, more 
feedback resources are required in the relay-based system. On 
the other hand, the proposed scheme is less complex compared 
to the simultaneous transmission-based relay-ARQ scheme. 
This is because the simultaneous transmission schemes, e.g., 
El-ED, are based on the assumptions that either different 
frequency bands are used by the source and the relay or the 
STC-based techniques are used while the source and the relay 


— RTD, optimal power allocation 
—INR, optimal power allocation 



Outage probability for 
single-user and relay—assisted 
networks, V‘^=V'^=0.5, 


Input SNR 101ogjjj((|) ) (dB) 


Figure 8. Outage probability for single-user and relay-assisted networks, 
sum power constraint (scenario 1), M = 1, initial transmission rate R = 0.5. 
For both single-user and relay-assisted channels, we set A®** = 1. For relay- 
assisted channel, A'*'^ = A"^"* = 0.5. 


are perfectly synchronized. With our data transmission model, 
none of these assumptions are required. 

On the optimal power terms: Fig. 10 shows an example 
of optimal power terms minimizing the outage probability. 
Different behaviors are observed at low and high SNRs, which 
can be interpreted as follows. 

At low SNRs, the power resources are limited and the ARQ 
scheme is conservative. In this case, Pf is set to be high, such 
that either the relay or the destination can decode the data at 
the end of round 1 and, as a result, P| is low. Then, if the 
relay decodes the message, high power is given to the relay 
to exploit the good properties of the relay-destination link and 
increase the probability of successful decoding. 

At high SNRs, i.e., when the sum power constraint is 
relaxed, there is enough power to gamble. Here, P* is set 
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Figure 9. Outage probability for different data transmission schemes, INR 
ARQ protocol, R = 0.5, M = 1, = 1, A"" = A""* = 0.5. 



Input SNR 101ogj|j((|)““^) (dB) 

Figure 10. Optimal power terms minimizing the outage probability of the 
relay-RTD network, /? = 1, M = 1, A*"* = 1, A'^'' = A''"* = 0.5. 


to be low (but high enough such that the relay can decode 
the message with high probability). If the source-destination 
channel is bad, the data can not be decoded by the destination 
and is retransmitted with higher powers. On the other hand, 
if the source-destination channel experiences good conditions, 
this gambling brings high return. Moreover, we have 
because 1) the relay-destination link experiences better average 
characteristics than the source-destination link and the relay 
requires less power than the source to guarantee a given outage 
probability. Also, 2) with high SNR, it is more likely that the 
data is retransmitted through the relay than the source. Thus, 
to keep the average sum power limited, less power is given to 
the relay than the source in the second round. 

On the effect of imperfect feedback signals: In wireless 
networks, the feedback signals reach the transmitter through a 
communication link experiencing different levels of noise and 
fading. Hence, it is probable to receive erroneous feedback 


signals at the transmitter(s). Setting M = 1, Fig. 11 studies 
the effect of imperfect feedback channels on the performance 
of the single-user and relay communication setups using RTD. 
Here, we have used the fact that, following the same arguments 
as before and with M = 1, the total average power and the 
outage probability of the relay-RTD network are respectively 
found as 

77/ 'Z'totall 

^total _ f 

- £:{ 7 -total} ’ 

^{ 5 total| ^ ps ^ ^ p ps) 

+ (1 - u;i)(l - p^PI) ((1 - /?i)(l - 7‘‘)(1 - + /3l7V‘') 

+ - ff) + P^CUi/3i), 

p{ 7 -total} ^ 

1 + (1 - (1 - cui)p-) ((1 - /?i)(l - 

+ (/3i + (1 - /3i)cui) - /3i + , 

(36) 

Pr(Outage) = (1 — 

+ (1 - - /‘>)(cci(l - 

-f (1 - wi)(l - p'^)a -f (1 - ^i)(l - a;i)p“), 

p = Pr(log(l+5^‘*P^+5^‘‘P^) <P), 
a = Pr(log(l+/‘>(Pi^ + P2*)) <P), 

= Pr(log(l + ff^Pl + Z)<R),Z= 

(37) 

if the feedback links are noisy. Here, p’^‘* and p*"^ denote 
the feedback bit error probability in the destination-source, 
the destination-relay and the relay-source feedback links, 
respectively. Moreover, the random variable Z = + 

in dJTl i comes from the fact that, with an erroneous 
feedback, the relay and the source may retransmit the data 
simultaneously, for instance, if the ACK bits sent by the 
destination and the relay are not correctly received by the 
source and relay at the end of the first round. Note that setting 
psd _ prd _ psi — Q jjj (|36l i- (|37l i. the results are changed to 
the ones in ([TtI i with M = 1. 

As shown in Fig. 11, the effect of imperfect feedback 
signals on the outage probability of the single-user and relay- 
assisted ARQ schemes is negligible, if the feedback bit error 
probabilities are in the practical range of interest (The same 
point is valid for the throughput, although not demonstrated 
in the figure.). However, the erroneous feedback signal affects 
the system performance at high feedback error probabilities. 

On temporal variations of the fading coefficients: In Figs. 
12-13, we demonstrate the outage probability and the through¬ 
put of the relay-ARQ setup in fast-fading conditions and the 
results are compared with the ones achieved in a quasi-static 
channel. Compared to the case with a quasi-static channel, 
better data transmission efficiency is observed in the fast¬ 
fading model. This is because with fast-fading more time 
diversity is exploited by the ARQ protocols. Also, compared to 
the quasi-static model, the effect of optimal power allocation 
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Figure 11. On the effect of noisy feedback channel. Optimal power allocation, 
RTD HARQ protocol, R = 0.5, M = 1, = 1, SNR = lOdB. For the 

relay-assisted setup, we have = A^^^ = 0.5. 



Input SNR 101og^^(f“'^) (dB) 


Figure 13. The throughput of the relay-ARQ setup in different fading 
conditions, INR protocol, R = 0.5, M = 1, A'^'^ = 1, A^*" = A*"^ = 0.5. 
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-^RTD, uniform power allocation 
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— RTD, optimal power allocation 
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Figure 12. The outage probability of the relay-ARQ setup in different fading 
conditions, R = 0.5, M = 1, A®'* = 1, A“ = A'"'* = 0.5. 


on improving the outage probability and the throughput of the 
relay-INR (resp. relay-RTD) setup increases (resp. decreases) 
when the channel is fast-fading. Moreover, the simulations 
show that, with fast-fading, the difference between the perfor¬ 
mance of the single-user and relay-assisted networks decreases 
slightly. However, there is still remarkable difference between 
the performances of these two communication setups and 
implementation of the relay leads to significant performance 
improvement. Finally, it is worth noting that we have tested 
the results of the fast-fading condition in many different cases, 
but since the results follow the same trend as the ones in the 
quasi-static model, we have chosen not to include them in the 
paper. 

Finally, it is interesting to note that, in harmony with El, 
EH, ca, in all cases better system performance is achieved 
by the INR ARQ protocol, compared to the RTD. 


V. Performance analysis in spatially-correlated 
LADING conditions 


This section studies the system performance in spatially- 
correlated conditions. As demonstrated, the discussions of 
the previous sections are helpful for performance analysis in 
spatially-correlated relay-ARQ setups. 

For simplicity, we concentrate on the scenario where the 
relay is close to the destination, while the source is far 
from them. As an example of such a case, we can consider 
the moving-relay systems, e.g., jn, m- Also, the same 
discussions are valid for the scenario with the source and the 
relay close to each other. 

Assuming the relay to be close to the destination, the 
relation between two fading random variables and is 
modeled by 


+ x/l - ^ 




(38) 


Here, c is a complex Gaussian variable CJ\f{0, j) uncorrelated 
with Also, (5 is a known correlation factor which demon¬ 
strates the dependency of the fading realizations. Moreover, A 
is the fading parameter of the source-relay and the source- 
destination links. This is a well-established model for the 
correlated Rayleigh-fading channels Il62l - ll64l . Under this 
model, the joint and marginal pdfs are found as 


fg^n«ix,y) = 


A2 


-A, 


1-S^ 


^Boi 


2X5y/Xy 

1-5^ 


), (39) 


fg^<{x) = Ae fg^i{x) = Ae (40) 

respectively, where Bq is the zeroth-order modified Bessel 
function of the first kind 1^ . Also, we have fgri{x) = 
^rdg-A which is independent of the two other variables. 

To study the effect of the spatial correlation, the key points 
are that 1) the two first steps of finding the closed-form 
expressions ©-(l^l are independent of the fading model. 2) 
The only terms which are affected by the fading model are 
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the probabilities Pr(Am), Pr(i?n,m), Pr(<S'm)- Moreover, these 
probabilities can be represented as tt = Pr(g*'' € n 

^sd g j^sd.i^ ^sd,2j PI ^rd g jyd,i^ ^rd,2j^ pj-Qper Selection of 

the boundaries Thus, to take 

the spatial correlation into account, it is only required to 
calculate tt as 




= &rd,l 


fg,i[x)y{x)Ax, 


y{x) = Pr e n e 


9^^ = X 


(41) 


where y{x) is determined based on the following procedure 



= Y (u, V, w, z). 

(42) 

Here, (a) is obtained by defining q = = y/WqS 

and using variable transform g = y^2y/q. Then, (6) is 
directly obtained from the definition of the Marcum Q-function 
M{x, y) = te ^ So(a:f)df. Finally, (c) is based on the 

fact that M.{x,y) = 1 + ^/‘^Boixy) — M{y,x) and 

[d) is derived by using variable transform t = y/x, partial 
integration and some calculations. 

As an example of dril l, consider the relay-RTD scheme in 
spatially-correlated quasi-static channel condition where, using 
dH, (Eoli, dull, the terms am,/3m, 7m, are rephrased as 




Pr(g- G [ 


- 1 


- 1 


ps ’ Y^m-l p, 
L^i=l ^ i Z^i=i -f i 


] n 5 ^“ G [ 0 , 


-1 


E m p; 
2=1 


d) 


git — 1 — 1 — 1 

ps , Y^m ps , _m —1 ps 

Z^i=l Xi 2^i-i xl 


(43) 


pR _l — 1 — 1 

R = Y(— _ —_ 0 —_) 

Rm ^ Vy^m ps, „m-l ps ’ ’ ^m-l p, 7 


(44) 


git _ 1 ^ - 1 


7m Y{Q, m ps’^’y^M ps)’ 
2^2=1 


(45) 


^T=l 


rd„-A'“x 


A“e 


Y 






y^i ps 


pr 

^ X Z^i=j+1 Xi 

y^i ps 
Z^i=l Xi 


- 1 


cR 


- 1 


E 3 ps' y^i-1 ps 
i=i Xi Z^i=i Xi 


dx. 


(46) 


Here, (l46l l is obtained numerically with Y{u, v, w, z) given 
in (l42li . Note that (|46] | is one-dimensional integration with 
known integration boundaries. Thus, m is calculated easily. 
Also, the other probability terms, e.g., i?„, are obtained with 
the same procedure as in (l46l l. Then, having the probabilities, 
the system performance is studied with the same procedure as 
before. 

Considering (l38]l . Figs. 14-15 show the relative coverage 
region and the outage probability of the relay-ARQ setup 
for different correlation conditions. The results indicate that 
in the practical range of correlation conditions the fading 
dependencies do not affect the system performance consid¬ 
erably, in the sense that the coverage region and the outage 
probability changes are negligible at low (5’s (The same point 
is valid for the throughput although not demonstrated in the 
figures). On the other hand, the data transmission efficiency of 
the relay-assisted network is considerably reduced at highly- 
correlated conditions, i.e., with <5 ~ 1. Specifically, the relay 
network is mapped to the source-destination single-user setup 
if (5 = 1 (please see (l38T l). The coverage region increases 
with the number of retransmissions substantially (Fig. 14). 
Moreover, although not demonstrated in the figure, the results 
indicate that the sensitivity of the coverage region to the fading 
parameter A“^ (resp. correlation factor S) decreases with <5 
(resp. input SNR). 

Finally, note that the performance gain of the relay-ARQ is 
at the cost of coordination overhead between the source and 
the relay. Particularly, the source data retransmission is decided 
based on the feedback signals from the destination and the 
relay, as opposed to non-relay setups with feedback only from 
the destination. Moreover, in harmony with every cooperative 
system, the mismatches between the source and the relay and 
the feedback delay may affect the data transmission efficiency 
of the relay-ARQ protocols. However, as shown in Fig. 11, 
the system performance is (almost) insensitive to the imperfect 
feedback signals for the practical range of the feedback error 
probabilities. 

VI. Conclusion 


In this paper, we studied the performance of the relay-ARQ 
networks using adaptive power allocation. The throughput and 
the outage probability were analyzed with different sum and 
individual power constraints on the source and the relay, under 
independent and spatially-correlated fading conditions. The 
results show that considerable outage probability reduction 
and throughput/coverage region increment are achieved by 
optimal power allocation between the source and the relay. 
The performance of the relay-ARQ network is not sensitive 
to spatial correlation, within the practical range. Also, the 
effect of imperfect feedback signals on the data transmission 
efficiency of the relay-ARQ scheme is negligible, if the feed¬ 
back bit error probability is in the practical range of interest. 
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Figure 14. Relative coverage region, compared to the single-user setup 
without ARQ, vs the correlation factor S in Input SNR=5dB, R = 0.5, 

e = 10-2, = 0.2. 



Figure 15. Outage probability vs the coiTelation factor S in (38). A = 1, 
= 0.2. RTD scheme with R = 0.5, M = 1. 


Finally, substantial performance improvement is achieved by 
variable-length coding and increasing the number of ARQ- 
based retransmissions. 


m—1 






h\og{l+g^^P:) ^ „ 

Y^m-1 , ^ ^(m-1) 


nE 


k\og{l+g^Pl) 


> R{m) n E 


U\og{l + g^^Pl) 


E m 1 — -^^{m) ' ' / ^ T^m j 

= Pr(V5-,5^V(m-i)(5“)<^(m-i) 

n > Rim) n 


^ Rim)) 


(47) 


for the INR. Here, we have defined 


^Im) = arg {[/(,„) (x) = R)m)}: 


(48) 


arg {[/^^ (x) (49) 

X 

Also, with the same arguments, the terms f3rn, 1 m and ojj in 
(ITSl) and (I 2 TI 1 are obtained by 

\ sd d \ sd d \ sr r 

/3^ = (e“^ - e"^ _ g-^ 

(50) 


= {1- (51) 


ujj = , (52) 

when the INR ARQ is implemented. Note that U)^^(O) = 0 
and (0) = 0. Moreover, as i?(j) < it is 

straightforward to show that and *l22l> 

and ( I 25 ] ), are increasing functions of x. Therefore, the solutions 
of (l48])-(l49b. i.e., and x'^^ will be unique for any 
given values of R(^rn), = 1, ■. ■, m. Unfortunately, to the 

best of authors’ knowledge, there is no closed-form solution 
for x^^^j and x^,^ . However, because of the uniqueness 

property, the solutions of (l48ll-(l49ll are easily obtained via nu¬ 
merical methods, such as the “fsolve” function of MATLAB. 

Finally, using the INR, the probabilities Oj^ra and pn should 
be calculated based on, e.g.. 


VH. Appendix 

A. Deriving the probability terms for the INR protocol in 
quasi-static channels 

Considering the INR, the probabilities 
Pr(Am), Pr(i?„ m), Pi'(<S'm) are determined with the 
same procedure as for the RTD protocol, i.e., (fT8T l- (l22l i. while 
the terms am, Pm, 1m, ^j, (^j,m and are recalculated 
based on the maximum decodable rate functions of the INR 
protocol. For instance, using (I23t- (l24l i. the probability in 
(HI is rephrased as 


/ / e-(^“‘’"+^''^)dxdy, (53) 

which is a two-dimensional numerical integration and, con¬ 
sequently, difficult to find. This is because the boundaries of 
the two-dimensional integrals can not be expressed in closed- 
form. To tackle the problem, two different methods can be 
used. The first one is the bounds introduced in Theorems 1-2. 
The second method is to use log(l -f x) ~ x for low SNRs 
which leads to the following approximation 
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p„ = Pr ( log(l + 

+ E ((lE - 7^)108(1+/“^^)) <1 


Jo 


1 y - ,pM ,_i _1 \ ps 

Z^i=n+1 V/{(j) R^i-i)'^i 




= 1 - 


1 -Er=i «1J^-1I7:^)1=|)» 

^rd^ ^(t) ^(i —1) 

^ V-M-/ I 1 5r,S ^ 

1 - e TTcTTiy) i I 






)^1) 




(i) «(i-l) 




1 - 






with Im representing the mxm identi^ matrix and C = \ci^k\ 
denoting the diagonal matrix given b}|j 


^i.k — 


^srp_s if i = k,i = 1,... ,m, 


0 


if i y k. 


(57) 


Using the Minkowski’s inequality 11661 Theorem 7.8.8] in 
leads to 


T' = det(I™+C) > (l + det(C)-)™. (58) 

Thus, from det(C) = (p*”)™ have dt > (1 + 


VTiZ^r and 

Pr(f]log(l + p^'/70<i?) < 




Pr 1 + 


^ 1 
Cm — 1 




(59) 


as stated in the theorem. 

For the second inequality of the theorem, i.e., (|29]l, the same 
arguments as in (l55]) -(l59b are used to write 

( j m 

J2 iog(i + g^^PD + E iog(i + ^ ^ 

i—1 i—j-\-l 

= Pr(T < e^), 

j m 

T = n(l + 5^'^') n (l + 3^‘’^D=det(I™ + D), 


i=j+l 


(60) 


( 54 ) where D = [d^^n] is the mxm diagonal matrix defined by 


The same approach can be used to find 6j^rn for low SNRs. 
Finally, having Pr(Am), Pr(i7„_m), Pr(>S'm) as functions of 
R{m), Pm, Py, m = 1,..., M + 1, the system performance 
can be studied with the same procedure as before (please see 
subsection IV.D). 

B. Proof of Theorem 1 

Considering (fTSTl-dZTli. it can be easily shown that the perfor¬ 
mance of the fixed-length INR scheme is a decreasing function 
of the probability terms Pr(X]™ i log(l + g^’^Pl) < R) and 

Pr(ELi iog(i + (fPl) + TZj+i iog(i + g^^Pl) < R)- In 

other words, the system performance is underestimated if the 

• j j 1 1 rri', fixed-length j rrd, fixed-length 

maximum decodable rates and are 

replaced by their corresponding lower bounds. From (|26] |. we 
can write 

Pr < ^) = Pr (^E ^ 


if fc = n,/c = 1,.. .,7, 

dk,n={ 9'‘^Pi if fc = n, fc = 7 -f 1,... ,m, 

0 if fc 7^ n. 

In this way, we reuse the Minkowski’s inequality to write 
T > (1-f det(D)^)’" 


(61) 


which, from the definition of Vj^mip) in (l29T l, leads to 
Pr(T < e^) < Pr ( < 


em-1 


VnLiTin™i+iPi 
= /o“ /g*‘'(a^)Pr(/‘* < -^)dx 

X 


sd r°° g-A“a:-A'‘’a:i-"*S(J^ 


= Pr(T' < e^). 


(55) 


= 1 ^ — 


R 

em -1 


Here, ^ is defined as 

m 

^ = n + 9"PD = det(I™ + C) (56) 




( 62 ) 


^The matrices are presented by capital bold letters. 
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C. Proof of Theorem 2 

To prove the theorem, the following modifications are 
applied into the arguments of Theorem 1. Considering ( l56l l, 
we rewrite Ti as 


Thus, the probability ^ lower bounded 


by 


Pr(T < e«) > Pr (^(l + ' 

™ r\/ i/r-1 / ! -^-\ 

'S = l[{l + g^^P:)=det{lm + G^^P^). (63) =j /^^x) Pr < V —^(1 + - ij dx 


Here, 0*“^ = p^'Im and P'’ = where C is given in (l57l) . 
By setting B = Im,A = X = P'^, denoting the conjugate 
transpose of the matrix X by X* and because the matrices 
and P® are Hermitian, we use 

det (AA* + BB*) det (I^ + X*X) > (det (B + AX)) ^ (64) 
lICTl Theorem 3.4] to write 

4- < v'det((G^^)2 +1™) det(I™ + (P'^)^) 


\ 


Z=1 




> Pr 




- ^9”(t 


(l + (p-)2)-J](l + (i^02)<ei 


- 1 ) 




-A" 


= 1 — e 




-1 


/d if fc = n,fc = 1,... ,j, 


9Z={ ifk = n,k = j + l,...,^ 

0 if fc n. 


if fc = n,fc = 1,... ,j, 

Pfc'n = ■^ P'i if fc = n,fc = j + 1,.. .,m, 
0 if k n. 

Then, we reuse ( l64t to write 


T = n n (1 + = det(B + G^‘''‘V'‘) 

2=1 


3 / \ PP^-3 

<\l[i + {9ZZ i + C, 


c=n n • ( 67 ) 

2=1 ^ ^ «^i+i ^ ^ 








^2R 




-A" 


■.j/f(l+x2)i-m -1 

( 68 ) 


(l + (gsr)2)™n(l + (i^02), (65) 


ULia + ip:r)UT=j+ii^ + iPD^) 

if the transmission powers are so low (or the initial transmis¬ 
sion rate R is so high) that r > 1. 

D. Proof of Theorem 3 

To obtain (l33l l. and the decoding probabilities of the RTD, 
let us define Zj^m = YZLi with 


Therefore, a lower bound of the probability 
Pr (fffJiLi log(l + 9^'^Pi) ^ P) is obtained by 

Pr(f]log(l+g^^i^^)<i?) 


f gfP- ifi = l,...,j, 

— A rd nr ■ ■ i i 

I 9iPi if 1 = ] + 

We have 

j m 

Pr(log(l + ^ Plff + ^ Plcft) < R) 


(69) 


i=l 


i=j+l 


.e«-l 

= PriZj^rn < - 1) = / „(x)dx, 

Jo 


(70) 


( 66 ) 


For (inil, i-e., the second inequality of the theorem, we redefine 
B = I„, A = G"‘’’''^X = F’'' such that and 

P"^=b^:n]with 


where fz ^ is the pdf of the random variable Zj^m ■ As the 
pdf of the sum of independent random variables is obtained 
by the convolution of their pdfs, we use ( 1691 ) and the inverse 
Laplace transform £“t| | to write 

fZj^m i^) — P { , P!s-,-nm , PJsx } 


■nLi(i + ^)nr=,+i(i+^)‘ 


of 


— 1+V + 


•? ps^sd 


2=1 


sd 

a, = 


Xsd 


E 

1 


PJaf -4a 




af = 




ni=i(1 x^p' ) T\.k=j+i,k=^i (f pr) 


(71) 

Here, (a) is based on (l6^ for independent Rayleigh-fading 
variables g^^ g'f^ i = 1,..., m, and 
1 




if i = 


i+- 


pT^ if i = j + 1,..., m. 


(72) 


Also, (b) is obtained by partial fraction expansion of V{s) = 

— - f PI PI f Plf + fc, with 

nLi(i+^)nr=.+i(i+^) * ^ ^ ^ ^ 
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and a’f, i = 1,..., m, representing the fraction expansion 
coefficients. Replacing ( fTlI i into (fTOl i leads to (l33l l, as stated in 
the theorem. Note that (fTlI i is based on the assumption that the 
function D{s) consists of m first-order poles, which is the case 
with optimal power allocation. Straightforward modifications 
should be applied in (l33T l and (fTlI i. if D{s) has poles of order 
n > 1. 

Equation (I74l i. on the other hand, is a direct consequence 
of ll60l Corollary 2] which, due to space limit, is not repeated 
here. Finally, using the same procedure as in (l33l l and (17^ . we 
can find am, Pm, 1m, ^j, Pn and the probability terms 
of the RTD and INR protocols, respectively. 
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